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Abstract 

Ambient RF (Radio Frequency) energy harvesting technique has recently been proposed as a 
potential solution to provide proactive energy replenishment for wireless devices. This paper aims to 
analyze the performance of a battery-free wireless sensor powered by ambient RF energy harvesting 
using a stochastic geometry approach. Specifically, we consider the point-to-point uplink transmission of 
a wireless sensor in a stochastic geometry network, where ambient RF sources, such as mobile transmit 
devices, access points and base stations, are distributed as a Ginibre a-determinantal point process (DPP). 

The DPP is able to capture repulsion among points, and hence, it is more general than the Poisson point 
process (PPP). We analyze two common receiver architectures: separated receiver and time-switching 
architectures. For each architecture, we consider the scenarios with and without co-channel interference 
for information transmission. We derive the expectation of the RF energy harvesting rate in closed form 
and also compute its variance. Moreover, we perform a worst-case study which derives the upper bound 
of both power and transmission outage probabilities. Additionally, we provide guidelines on the setting 
of optimal time-switching coefficient in the case of the time-switching architecture. Numerical results 
verify the correctness of the analysis and show various tradeoffs between parameter setting. Lastly, we 
prove that the sensor is more efficient when the distribution of the ambient sources exhibits stronger 
repulsion. 
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I. Introduction 

Ambient RF energy harvesting techniques offer the capability of converting the received RF 
signals from environment into electricity [1], [2]. Therefore, it has recently emerged as an 
alternative method to operate low-power devices [3], [4], such as wireless sensors [5]. Ambient 
RF energy harvesting aims to capture and recycle the environmental energy such as broadcast 
TV, radio and cellular signals [6], which are essentially free and universally present, making 
this technique even more appealing. An experiment with ambient RF energy harvesting in [7] 
shows that 60//W is harvested from TV towers that are 4.1km away. It is also reported in [8] that 
109//W RF power can be harvested from daily routine in Tokyo. In [9], the authors measure the 
ambient RF power density from 680MHz to 3.5GHz and show that the average power density 
from 1GHz to 3.5GHz is of the order of 63// W/m 2 . Detected 6.3km away from Tokyo Tower, 
the RF-to-DC conversion efficiency is demonstrated to be about 16%, 30% and 41% when the 
input power is -15 dBm, -10 dBm and -5 dBm, respectively [10]. 

In this context, wireless devices powered by ambient RF energy are enabled for battery- 
free implementation, and a perpetual lifetime. For example, reference [11] demonstrates that an 
information rate of 1kbps can be achieved between two prototype devices powered by ambient 
RF signals, at distance of up to 2.5 feet and 1.5 feet for outdoors and indoors, respectively. 
Existing literature has also presented many implementations of battery-free devices powered by 
ambient energy from WiFi [12], GSM [13] and DTV bands [14] as well as ambient mobile 
electronic devices [15]. 

A. Related Work 

Geometry approaches have been applied to analyze RF energy harvesting performance in 
cellular network [16], cognitive radio network [17], and relay network [18]-[20]. The authors in 
[16] investigate tradeoffs among transmit power and density of mobiles and wireless charging 
stations which are both distributed as a homogeneous Poisson Point Process (PPP). Energy 
harvesting relay network has been mostly analyzed. In [17], the authors study a cognitive radio 
network where primary and secondary networks are distributed as independent homogeneous 
PPPs. The secondary network is powered by the energy opportunistically harvested from nearby 
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transmitters in the primary network. Under the outage probability requirements for both coex¬ 
isting networks, the maximum throughput of the secondary network is analyzed. The study in 
[18] analyzes the impact of cooperative density and relay selection in a large-scale network with 
transmitter-receiver pairs distributed as a PPP. Reference [19] investigates a decode-and-forward 
relay network with multiple source-destination pairs. Under the assumption that the relay nodes 
are distributed as a PPP, the network outage probability has been characterized. The studies in 
[20] investigates network performance of a two-way network-coded cooperative network, where 
the source, destination and RF-powered relay nodes are modeled as three independent PPPs. 

Other than RF energy harvesting, stochastic geometry approaches have also been applied 
to address other types of energy harvesting systems. Reference [21] investigates the network 
coverage of a hexagonal cellular network, where the base stations are powered by renewable 
energy, and the mobiles are distributed as a PPP. The authors in [22] explore the network coverage 
in a relay-assisted cellular network modeled as a PPP. Each relay node adopt an energy harvesting 
module, the energy arrival process of which is assumed to be an independent and identical poisson 
process. In [23], the authors provides a fundamental characterization of the regimes under which 
a multiple-tier heterogeneous network with genetic energy harvesting modules fundamentally 
achieves the same performance as the ones with reliable energy sources. Different from above 
studies, our previous in [24] adopts a determinantal point process model to analyze the downlink 
transmission performance from an access point to a sensor powered by ambient RF energy. 

B. Motivations and Contributions 

As discussed in Section I-A, the prior literature mainly focuses on the performance analysis 
on RF-powered wireless devices using PPPs. We generalize this approach by considering a 
larger class of point processes, specifically the Ginibre a-determinantal point process. Table I 
compares the PPP and the Ginibre a- DPP with regards to a few key points summarized in 
the table. Therein, the column “simulation” refers to the ease (in terms of time, computational 
complexity and implementation difficulty) of the simulation of the point process. By mathematical 
tractability, we mean the possibility of obtaining closed mathematical formulas for the moments 
of the point process. By modulability, we mean the choice of available parameters. Lastly, by 
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TABLE I 

Comparison of common stochastic models. 


Model 

Simulation 

Mathematical tractability 

Modulability 

Data fitting 

PPP 

Very easy 

Closed forms 

Intensity 

Not very fitting 

Ginibre a-DPP 

Easy 

Closed forms 

Intensity, a 

Different degrees of repulsion 


data fitting we refer to the range of phenomena modeled by the point process. 

The Ginibre o-DPP offers many advantages in terms of modeling capability and ease of 
simulation [25] (here, —l<a<0isa parameter, and the PPP is a special case obtained in the 
limit a —> 0). One advantage of the Ginibre o-DPP over the PPP is that the Ginibre o-DPP can 
be used to model random phenomena where repulsion is observed. Mobile systems may exhibit 
some clustering and repulsion behaviors, such as in mobile sensor networks [26], mobile cellular 
networks [27] and mobile social networks [28]. Therefore, the Ginibre o-DPP is a suitable tool 
for the analysis of the impact of distribution patterns on network performance. Previous works 
in [29] and [30] have adopted Ginibre point processes to model the locations of base stations in 
wireless networks. 

Since the direct simulation of DPP models is computationally slow, in this paper we fo¬ 
cus on a worst-case scenario which is simpler to analyze. Namely, we perform a worst-case 
analysis of the point-to-point uplink transmission between an RF-powered sensor node and a 
data sink. The sensor node needs to harvest RF energy from ambient RF sources (e.g., cellular 
mobiles and access points), which are modeled using a Ginibre o-DPP which we have briefly 
described previously. The sensor, assumed to be battery-free, transmits to the data sink using 
instantaneously harvested RF energy. A power outage happens if the instantaneously harvested 
energy fails to meet the circuit power consumption of the sensor. Moreover, if the minimum 
transmission rate requirement cannot be fulfilled, a transmission outage occurs. For the study 
of the transmission outage probability, we consider the scenarios of out-of-band transmission 
and in-band transmission. In the former scenario, the sensor node transmits data on a frequency 
band different from that for RF energy harvesting (without co-channel interference). In the latter 
scenario, the sensor node transmits on the same frequency band of ambient RF energy sources 
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(with co-channel interference). We focus on analyzing the performance of two different receiver 
architectures: a separated architecture [1] and a co-located receiver architecture, called time¬ 
switching [31], whose details are in Section II-A. By modeling the RF sources with a Ginibre 
o-DPP. we analyze the impact of the distribution of ambient RF sources on the performance of 
the RF-powered sensor. Our main contributions are summarized below. 

• First, we derive the expectation of aggregated energy harvesting rate by the sensor. We 
also obtain the expression of the variance of the energy harvesting rate. The closed-form 
expressions are verified by the numerical results. Since obtaining numerical results is very 
time-consuming, our closed forms are very useful in practice. 

• Next, we investigate the power outage probability, i.e. the probability that the sensor node 
becomes inactive due to lack of sufficient energy supply. We give an upper-bound of the 
power outage probability in closed form, and we interpret this upper-bound as a worst-case 
scenario. We confirm that the theoretical power outage probability and its estimation by 
simulation are consistent. We also compare the performance of separated and time-switching 
receiver architectures by simulation. 

• We further study the transmission outage probability, i.e. the probability that the sensor 
fails to fulfill its transmission rate requirement, because of insufficient transmit power. The 
scenarios of out-of-band transmission and in-band transmission are both considered. We 
derive an upper-bound of the transmission outage probability and again interpret is as a 
worst-case scenario. Inspired by the observations in the numerical results, we derive the 
optimal value of the time-switching coefficient r. The expression for the optimal choice of 
r is given and verified by numerical computation. Lastly, we derive a lower bound of the 
transmission rate. 

Note that this paper is an extension of [32], wherein we present partial of the results with the 
separated architecture. 

The remainder of this paper is organized as follows. Section II introduces the system model, 
the DPP geometry model of ambient RF sources and performance metrics. Section III estimates 
the performance metrics of the sensor for both Ginibre o-DPP and PPP modeling of ambient 
RF sources. Lastly, our conclusion can be found in Section IV. 
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Fig. 1. A network model of ambient RF energy harvesting. 


Notations : Throughout the paper, we use E[X] to denote the probabilistic expectation of a 
random variable X, and P(A) to denote the probability of an event A. 

II. System Model 

A. Network Model 

We consider a network comprising a number of ambient RF energy sources, i.e., wireless 
information transmitters, in which a sensor node is powered solely by the energy harvested 
from these energy sources. Figure 1 shows the considered network model, where the sensor 
node harvests RF energy emitted from the ambient sources and utilizes the harvested energy 
to perform uplink data transmission to the data sink. We model the distribution of ambient RF 
energy sources as a Ginibre ck-DPP, which will be specified in detail in Section II-B. The transmit 
power of the ambient RF sources are assumed to be identical. Without any loss of generality, 
the sensor is considered to lie at the origin. Furthermore, we assume that the sensor node is 
battery-free. In particular, the sensor utilizes the instantaneously harvested RF energy to supply 
its operations. We study two different receiver architectures: separated receiver architecture, and 
time-switching, which either enables the sensor to perform data transmission and RF energy 
harvesting simultaneously or separately. 

• Separated receiver architecture: As shown in Fig. 2a, this architecture equips the energy 
harvester and the information transmitter with separated antennas so that they can function 
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(a) Separated Receiver Architecture (b) Time Switching Achitecture 

Fig. 2. Separated receiver and time switching architectures. 


independently and concurrently. The instantaneously harvested energy is first used to operate 
the sensor circuit and then the surplus energy is provided for information transmission. This 
architecture can maximize the utilization of energy harvesting devices, but it is generally 
larger in size compared to the time-switching architecture. 

• Time-Switching Architecture: As shown in Fig. 2b, the time-switching architecture, is 
equipped with a single antenna. By adopting a switcher, this architecture allows either the 
energy harvester or the information transmitter attached to the antenna at a time. The time¬ 
switching architecture works on a time-slot basis. In each time slot, the energy harvester 
first uses r (0 < r < 1) portion of a time slot to harvest RF energy. The capacitor reserves 
the surplus of the harvested energy after being used to power the sensor circuit. Next, 
during the rest of 1 — r time, the information transmitter utilizes the surplus energy from 
the capacitor to transmit information. 

The RF energy harvesting rate of the sensor node from the RF energy source k in a free-space 
channel P $ can be obtained based on the Friis equation [33] 1 as follows: 


Ph = QPPs 


G.GuX 2 


( 1 ) 


(4vr4) 2 ’ 

where /3 is the RF-to-DC power conversion efficiency of the sensor node, and g is an efficiency 
factor which depends on the specific architecture. For a given RF energy source, Ps is its transmit 
power, G g is its transmit antenna gain, A is the wavelength at which it emits. As the focus of 


'Other RF signal propagation models can also be used without loss of generality in the analysis of this paper. 



























this paper is to analyze the impact of the locations of ambient RF sources to the performance 


of the sensor node, we intentionally make other parameters, e.g., P s , G s , and A to be constants 
for ease of presentation and analysis. Nevertheless, the proposed analytical framework can also 
be extended to the case when these parameters vary. d k is the distance of an RF energy source 
k to the receiver antenna of the sensor node. Gu is the receive antenna gain of the sensor node. 
Let X/,. e M 2 be the coordinates of the RF energy source k (recall that the sensor node lies at 
the origin). The distance is modeled as d k = e + 1111 , where e is a fixed (small) parameter 
which ensures that the associated harvested RF power is finite in expectation. Physically, e is 
the closest distance that the RF energy sources can be to the sensor node. 

The aggregated RF energy harvesting rate by the sensor node from the ambient RF sources 
can be computed as 



( 2 ) 


where /C is a random set consisting of all RF energy sources. We assume that /C is a point 
process [34]. 

The sensor consumes a base circuit power, denoted by Pq. Note that this circuit power 
consumption also accounts for the energy loss due to various factors such as capacity leakage. 
Following practical models [35], the circuit power consumption of the sensor is assumed to 
be fixed. We assume that, other than circuit power consumption, there is no power loss during 
transfer from energy harvester to information transmitter for both architectures. For the separated 
receiver architecture, the transmit power is given by P T = [P H — Pc] 4 , where [.x] + = max(0, x ) 
and Pq is a constant. For the time-switching architecture, all the harvested energy can be used in 
the data transmission phase to maximize transmission rate. Thus the transmit power is dependent 
on the transmission time, and is given by Pj = \I\\ — Pc] + /(I — r). Then, the general form of 
maximum transmission rate of the sensor node is given as follows 2 : 



(3) 


Note that state-of-the-art wireless information receivers are not yet able to achieve this rate upper bound due to additional 


processing noise such as the RF band to baseband conversion noise. 




where W is the transmission bandwidth, and 0 < // < 1 is an efficiency factor depending on the 
specific architecture, a 2 is a nonnegative constant which represents the power of additive white 
Gaussian noise (AWGN). The term CYlkeic ^h corresponds to the interference, and the specific 
value of £ e {0,1} depends on whether we consider an out-of-band or in-band transmission 
scenario. h 0 denotes the channel gain between the transmit antenna of the sensor node and 
the receive antenna of data sink. The separated receiver architecture corresponds to g = 1 and 
r] = 1. The time-switching architecture corresponds to the g = r and p = 1 — r, where r is the 
time-switching parameter. In both of these cases, £ = 1 corresponds to an in-band transmission 
scenario, while £ = 0 corresponds to an out-of-band transmission. 

B. Geometric DPP Modeling of Ambient RF Energy Sources 

As an extension of the Poisson setting, we model the locations of RF energy sources using a 
point process K. on an observation window OcR 2 such that 0 < \0\ < +oo (here \0\ denotes 
the Lebesgue measure of O). In other terms, K, is an almost surely finite random collection 
of points inside O. We refer to [34] and [36] for the general theory of point processes. The 
correlation functions // n) of /C (if they exist), w.r.t. the Lebesgue measure on R 2 , verify 



E 


(4) 


for any family of mutually disjoint bounded subsets B 1 ,..., B n of E, n > 1. Heuristically, p <y> 
is the particle density, and (P l) (x \...., x n ) dx] ... dx n is the probability of finding a particle 
in the vicinity of each x,, i = 1 ,,n. The correlation functions are thus a generalization of 
the concept of probability density function to the framework of point processes. The correlation 
functions play an important role in the definition and interpretation of a general o-DPP. 

1) General a-determinantal point process: We let a = — 1/j for an integer j e N*, and we 
define a general a-DPP in the following. Let us introduce a map K : L 2 (M 2 ) i— * L 2 (M 2 ), where 
L 2 (M 2 ) is the space of square integrable functions on M 2 . We assume in the following that /C 
satisfies Condition A from [39]. The map K is called the kernel of the o-DPP It represents the 
interaction force between the different points of the point process. A locally finite and simple 
point process on R 2 is called an o-DPP if its correlation functions w.r.t. the Lebesgue measure 
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on R 2 exist and satisfy 

P {n \x i, ...,x n ) = det a (K(xi, Xj))i<ij< n , (5) 

for any n > 1 and x\, ..., x n e M 2 , and where the a-determinant of a matrix M = (Mjj)i<y< n 
is defined as 

n 

det Q M = II (0 > (6) 

crSSn 1=1 

where S n stands for the n-th symmetric group and u(a) is the number of cycles in the permutation 

a e S n . 

Let us now give some basic properties of the a-DPP to emphasize the role played by the 
kernel K. We start by a proposition exhibiting the repulsion properties of the a-DPP. Its proof 
follows from the formula defining the correlation functions (4). 


Proposition 1 (Repulsion of the a-DPP). The covariance of an a-DPP of kernel K is given by 

Cov(1C(A), )C(B)) = a f \K(x, y)\ 2 dxdy, 

JAxB 

where fC(A) and JC(B)) denote the random number of point process points located within the 
disjoint bounded sets A, B C R 2 . 

Since a < 0, 1C (A) and K,(B) are negatively correlated and the associated a-DPP is known to 
be locally Gibbsian, see, e.g., [38], therefore it is a type of repulsive point process. Additionally, 
the a-DPP exhibits more repulsion when a is close to — 1. As a —> 0, 1C (A) and JC(B) tend 
not to be correlated, and in fact it can be shown that the corresponding point process converges 
weakly to the PPP, c.f. [39]. 

Next, we recall from [40] the following proposition which gives the hole probabilities of the 
a-DPP. Proposition 2 allows us to compute the quantities known as hole probabilities. 

Proposition 2 (Hole probability of the a-DPP). For every bounded set B cR 2 we have 

P(/Cnfi = 0) = Det(Id + aAB)- 1/a , (7) 

where Kb is the operator restriction of K to the space L 2 (B) of square integrable functions 
on B with respect to the Lebesgue measure. Here, Id is the identity operator on L 2 (B) and for 
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any trace class integral operator K, Det (Id + aK) is the Fredholm determinant of Id + oK 
defined in [37]. 

2) The Ginibre point process: In the rest of the paper, we focus on the Ginibre o-DPP. which 
is a particular o-DPP well-suited for applications. The Ginibre process is a type of o-DPP that 
is invariant with respect to rotations. Therefore, it will be fruitful for computational convenience 
to restrict our attention to the choice of observation window O = B((). IT), defined as a disc 
centered around 0 and of radius R > 0. 

The Ginibre process is defined by the so-called Ginibre kernel given by 

K(x,y) = / , e ^ e -“ (|l|2+l!,|2) , x,yeO = B(0,R), (8) 

where p > 0 is a fixed parameter called density of the point process. This kernel is that of the 
usual Ginibre process defined, e.g., in [25], to which we have applied a homothety of parameter 
y/Wp > 0: x i y x/{y/Wp). Next we recall a few features of the Ginibre process. 

• The Ginibre kernel K defined in (8) satisfies Condition A from [39], and is thus a type of 
a -DPP. 

• The intensity function of the Ginibre process is given by 

p (1) (x) = K(x,x) = p, (9) 

c.f. [39]. This means that the average number of points in a bounded set B C B(0,R) 
is p\B\. Note that the intensity function of a homogeneous PPP is also constant, so p is 
interpreted as the intensity of the corresponding PPP. 

• The Ginibre o-DPP is stationary and isotropic, in the sense that its distribution is invariant 
with respect to translations and rotations, c.f. [25]. Hence, the Ginibre point process models 
a situation where the RF energy sources are distributed homogeneously in M 2 . 

We write K, ~ Gin(a, p) when /C is an o-DPP with Ginibre kernel defined in (8) and density p. 
The spectral theorem for Hermitian and compact operators yields the decomposition K(x. y) = 

where (<£>*)»>0 is a basis of eigenvectors of L 2 (0), and (Aj);> 0 are the 
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corresponding eigenvalues. In, e.g., [25], it is shown that the eigenvalues of the Ginibre point 
process on O = 13(0, R) are given by 

T(n + 1 , irpR 2 ) 


^ M 


n\ 


n E N, 


where 


r(z, a) — / e t z 1 df, z6 C, a > 0, 


( 10 ) 


( 11 ) 


is the lower incomplete Gamma function. On the other hand, the eigenvectors of K are given by 
ip n (z) = y=^=e (y/Wpz) n , for n E N and z E O. We refer to [25] for further mathematical 
details on the Ginibre point process. 

To illustrate how the parameter a affects the distribution of the DPP, in Fig. 3 we show some 
snapshots of the scattering of ambient RF energy sources in a disc of radius R = 10, when the 
RF source density is p = 0.3. It is seen that strong repulsion exists between the RF sources 
when a = —1. As a result, the RF sources tend to scatter evenly over the area. We can observe 
that the repulsion decreases very fast with the increase of a. When a = —0.5, some of the RF 
sources exhibit attraction by locating close to each other. Some grids of vacant area begin to 
emerge. The attraction keeps increasing as a approaches zero. When a = —0.03, the RF sources 
show clustering behavior, which is a feature of the PPP Consequently, there appear many grids 
of vacant area. Thus, depending on the distribution of RF sources that is observed, we shall 
choose a different value of the parameter a that appropriately models the situation at hand. 


C. Performance Metrics 

We define the performance metrics of the sensor node as the expectation of RF energy 
harvesting rate, the variance of the RF energy harvesting rate, power outage probability and 
transmission outage probability. Let us first introduce the mathematical quantities of interest. 

1) Theoretical values: The expectation of the RF energy harvesting rate is defined as Ep H = 

E[P h ]- 

The variance of RF energy harvesting rate is given by Vp H = E [(P H — E [Ph]) 2 ] • 

Power outage occurs when the sensor node becomes inactive due to lack of enough energy 
supply. The power outage probability is then defined as P po = P (P H < Pc) • 
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Fig. 3. Snapshots of the distribution of ambient RF energy sources (a) a = —1 (b) a = —0.5 (c) a = —0.03. 


Let m > 0 denote the minimum transmission rate requirement. If the sensor fails to achieve this 
requirement, a transmission outage occurs. The transmission outage probability can be defined 

as P to = P (C < m). 

2) Estimation by simulation: The different theoretical performance metrics introduced in 
Section II-C1 may in practice be estimated by Monte Carlo simulation of the underlying a- 
DPP. The simulation of a-DPPs when a = — 1/j, j e N, is done by using the Schmidt 
orthogonalization algorithm developed in full generality in [41], and specifically in [25] for 
the Ginibre point process. The simple generalization to a = — 1/j can be found in the recent 
survey [42], and additional details on DPP can be found in [43]. 

3) Upper bounds under a worst-case scenario: In practice, there are no closed forms for 
all the performance metrics appearing in Section II-C1. Additionally, estimation by simulation 
suffers from some drawbacks: 1) the time required to draw N samples of Pu can be rather long; 
2) the estimation may differ significantly from the theoretical result (if N is not sufficiently 
large); 3) it is difficult to modify the parameters retroactively, and a new set of simulations is 
then required. These drawbacks motivate the computation of upper bounds which is the object 
of the remainder of the paper. 

To that end, we introduce a simpler scenario, which is henceforth called worst-case scenario. 
In this worst-case scenario, the sensor node only receives energy from one RF source at a time. 
In this simpler case, there is energy (respectively transmission) outage if and only if the closest 
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(a) General Case: No Outage, Worst (b) General Case: No Outage, Worst (c) General Case: Outage, Worst Case: 
Case: No Outage Case: Outage (High RF Source Density) Outage (Low RF Source Density) 

Fig. 4. Example scenarios for ambient RF energy harvesting. 


RF source is further than some characteristic distance. Figure 4 illustrates the difference between 
the general-case scenario, and the worst-case scenario. 7 represents the maximum distance from 
the sensor node where a single RF energy source can still power the sensor node by itself. In 
Fig. 4a, as an RF energy source lies in the range of 7 , sufficient power is guaranteed from this 
single source. As a result, both the general-case scenario and worst-case scenario experience no 
outage. If there exists no RF energy source in the range of 7 , the sensor may still be powered if 
the sum of harvested energy from multiple RF energy sources is large enough. In this context, 
the outage largely depends on the density of RF energy sources. When the density is high, as 
in Fig. 4b, the sensor shall experience no outage in general-case scenarios. However, in the 
worst-case scenario, outage occurs. Contrarily, when the density is low, as in Fig. 4c the sensor 
does not harvest enough energy in both general-case and worse-case scenarios. 

As illustrated here, our worst-case scenario over-estimates the outage probability. Indeed in 
Fig. 4b, there is an outage in the worst-case scenario, although there is not in the general 
scenario. Therefore, the performance metrics in this worst-case scenario over-estimate the real 
performance metrics, this will be the focus of Section III in which we show that this worst-case 
scenario constitutes an upper-bound to the outage probabilities. 

III. Performance Analysis 

In this section we analyze the performance metrics defined in Section II-C when /C ~ Gin (a-, p) 
is the Ginibre a:-DPP with parameter a = — 1/j for j e N*, and density p > 0. 
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A. The Expectation and Variance of RF Energy Harvesting Rate 

First, we obtain the expectation of RF energy harvesting rate in the following theorem. 


Theorem 1. The expectation ofRF energy harvesting rate can be explicitly computed as follows 3 : 

e 


G^GuX 2 


E[P h ] = 2n Q pP s ^ A _l p 


■'e-tO 


(4vr ) 2 r \R + e 
pg/3PsGsGnX 


+ ln(P + e) — 1 — ln(e) 


8n 


2 In (j) . 


( 12 ) 

(13) 


Additionally, the variance of the RF energy harvesting rate can be computed as follows: 
nr 1 \ 2 \ 2 


V Ph = qPP s 


(4tt) 


27T P FT 


3R + e \ 2 

+ ap 


6e 2 6 (P + e) 3 

where recall that O = 13(0, R) is the obserx’ation window. 


r g— 7T p|l x ~y|| 2 

7 i II in 2 r + ii m 2 dxdy I ’ (14) 

oxo (e+ ||x||) (e+ ||y||) 


Before moving on to the proof, a few remarks are in order. First, we note that Theorem 1 
implies that at the level of expectations, the Ginibre a-DPP behaves like a homogeneous PPP 
and in particular, the expectation of RF energy harvesting rate is independent of the repulsion 
parameter a. Therefore, on average, the harvested energy is the same when a varies. However, 
it is straightforward from (14) that the variance of the RF energy harvesting rate is larger when 
the point process is closer to a PPP. Heuristically, there is a larger probability that there are no 
points close to the sensor when the RF sources are distributed as a PPP. Second, notice that the 
second term in (14) is in fact not a closed form. To the best of our knowledge, the second term 
cannot be explicitly calculated, but should instead be approximated numerically. 

Proof: We have 


by Campbell’s formula [34], where p^\x) = K(x,x) 
by (9). We thus find 


' + m 


= p is 


i > 

the intensity function of K. given 


E[P h ] = PPs 


G$Gu\ 2 


»R 


(4vr) s 


- 2 tt / p 


(e + r)‘- 


dr, 


> l. 


3 Here, we say that / ^ e _> 0 9 if //s 


e-MD 
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TABLE II 

Parameter Setting. 


Symbol 

/~ii s~ik 

p 

Ps 

W 

A/c 

Pc 

a 2 

Value 

1.5 

0.3 

1W 

lKHz 

0.167m 

15.8fj,W 

—90 dBm 


RF Energy Harvesting Rate versus Density of Ambient RF Energy Sources 
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Density of Ambient RF Energy Sources 


Fig. 5. RF Energy Harvesting Rate versus Density of Ambient RF Energy Sources. 


by polar change of variable, and the integral on the r.h.s. is computed explicitly as 

rR 


dr = 


R + e 


+ In (R + e) - 1 - ln(e) , 


Jo % + H 

which yields the result. 

For brevity, the proof of the expression of variance (14) is presented in Appendix I ■ 
Next, we examine the validity of the expressions of the expectation and variance of RF energy 
harvesting rate. All the network simulations in this paper are considered in the scenario of an 
LTE network with a typical 1800 MHz operating frequency. The corresponding wave length 
adopted is 0.167m. The channel gain of both transmit antenna and receive antenna are assumed 
to be 1.5. The RF-to-DC power conversion efficiency is assumed to be 30%. We consider the 
transmit power of ambient RF sources is 1 W. The circuit power consumption of the sensor is 
fixed to be -18 dBm (i.e., 15.8 fiW) as in [44]. The sensor is assumed to be allocated with a 
1 kHz bandwidth for data transmission. The AWGN is considered to be -90 dBm. The channel 
gain between the sensor and data sink is calculated as h 0 = 62.5 d~ A [45], where d is the distance 
between the sensor node and the data sink. The results with the separated and time-switching 
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receiver architectures are labeled as “SA” and “TS”, respectively. Note that the results for the 
PPP are identical to that of the a-DPP, when a — 0. Additionally, the performance of separated 
receiver architecture in terms of expectation and variance of RF energy harvesting rate and power 
outage probability is identical to the case when r = 1 for time-switching. 

As shown in Fig. 5, the numerical results, averaged over N = 5 x 10 5 of simulation runs, 
match the analytical expression (12) accurately over a wide range of density p, i.e., from 0.01 to 
1. As expected, the RF energy harvesting rate increases with the density of ambient RF energy 
sources. Under the same density, the RF energy harvesting rate is affected by e. We observe 
that when e = 0.001, larger RF energy harvesting rate is available at the sensor than that when 
e = 0.01. The straightforward reason is that, from (2), the smaller distance the RF sources locate 
near the sensor, the more aggregated RF energy harvesting rate can be achieved. Additionally, the 
use of the approximate expression (13) from (12) can be observed to increase with the density 
of ambient RF energy sources. However, the difference between (13) and (12) in percentage 
remains the same, i.e., 15 percent when e = 0.01. This difference is dependent on e, and is 
shown to diminish with the decrease of e. As shown in Fig. 5, when e = 0.001, the approximate 
expression approximates more closely to the analytical expression. 


B. Upper Bound of the Power Outage Probability 

In this section, we derive upper-bounds of power outage probability defined in Section II-C1. 
We interpret these upper-bounds in terms of a worst-case scenario, as specified in Section II-C3. 
We also point out that all the numerical estimations done in this section are performed under 
the worst-case assumption. 


Theorem 2. Let us define 


Then, the following bound holds: 


a_ A / q/3PsGsGh 

7 “ Ur V ' 


p — 

1 po 


H 


< Pc) < ( II ( 1 + Ot- 


\n> 0 


r(n + 1, Tip inf(f?, y) 2 ) 


n\ 


-1/a 


( 15 ) 


where F(z, a) is the lower incomplete Gamma function defined in (11). 
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Remark that the parameter e does not appear in the bound of Theorem 2 and rigorously, the 


inequality holds only for a sufficiently small value of e. In practice, we should therefore make 
sure that e is chosen small enough. 


Proof: To make the proof easier to follow, let us set /(x fe ) = o(5I\ , n- • f° r ^ e 


(47r(e+||x fe ||)) 2 ’ 


Then, 

P po = P /(x fc ) < P c J < P(Vfc e 1C, /(xfc) < P c ) = P(VA: e /C, ||x*|| > 7 - e) ( 16 ) 



P(/C D P( 0,7 — e) = 0), 


where we have chosen e such that 7 — e > 0. Thus by Proposition 2, we obtain 


Ppo < Det(Id + aiT B(0)7 _ e) ) 1/q . 


(17) 


Since in our case K is the Ginibre kernel, the eigenvalues of K are given by (10). By standard 
properties of the Fredholm determinant which can be found, e.g., in [37], we find, 



and the result follows by letting e go to zero on the r.h.s. of (17), since the associated function 


of e is continuous. 


Let us explain briefly the heuristics behind the bound in Theorem 2. Note that the only line 
which is not an equality in the proof of Theorem 2 is (16). The approximation made therein 
is that the harvested energy rate I \[ is provided by the RF source closest to the sensor node. 
Hence, the probability appearing on the r.h.s. of (15) is the probability corresponding to the 
worst-case scenario introduced in Section II-C3. Estimation of this upper-bound by simulation 
will therefore consist in assuming the worst-case scenario. Conversely, estimation of the l.h.s. 
of (15) will consist in assuming the general-case scenario. 

It should be noted that the eigenvalues appearing in the product of Theorem 2 are in decreasing 
order, and decrease exponentially when n>np inf(P, y) 2 , see [25] for details. Hence, the product 
which appears in Theorem 2 is well approximated by 



( 18 ) 
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Upper Bound of Power Outage Probability versus Density 
of Ambient RF Energy Sources 



Power Outage Probability versus p and x 



Fig. 6. Upper bound of power outage probability versus Fig. 7. Upper Bound of Power Outage Probability versus 
density of ambient RF energy source (separated receiver Density of Ambient RF Energy Sources and Time-Switching 
architecture). Coefficient r (time-switching architecture). 


where N 7 rpinf(i?, 7 ) 2 . We also note that 


dct 


In ( 


-1/a 


\n> 0 


1 (1 + a\ n ) ln(l + a\ n ) — aA n 

a 2 1 + a\ n 

n> o 


which means that the bound of Theorem 2 is smallest when a = — 1, i.e. when repulsion is 
maximal, and increases with a. 

As a corollary of Theorem 2, we find in the case of a PPP (which is obtained as the limit as 
a —> 0 in the theorem) the next corollary. 


Corollary 1 . Let 1C ~ Poiss(0,p) be a Poisson process on O = 13(0, R) with density p. Then, 
the following bound holds: 

P P o < e-^ inf ^ )2 , (19) 

where 7 is as defined in Theorem 2. 


In Fig. 6 , we show the variation of the upper bound of the power outage probability P po as 
a function the density of ambient RF energy sources p for both separated and time-switching 
architecture. It is observed that P po is a decreasing function of p. The numerical results verify 
that the analytical expressions for the upper bounds in (15) and (19) are very tight for different 
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values of a. We also observe that the smaller the a, the lower the upper bound of the outage 
probability. In other words, more repulsion among the locations of the RF sources results in 
lower outage probability of the sensor. This finding coincides with the previous analysis that the 
maximal repulsion results in best power outage performance. Moreover, we can observe that the 
influence of a depends on the density p. For both receiver architectures, the gap between the 
bounds for DPP (a = — 1) and PPP first increases with the density p and then decreases when 
the power outage probability becomes low. 

Figure 7 examines the impact of time-switching coefficient r on the upper bound of the outage 
probability. It is obvious that, regardless of density p, power outage probability is a monotonically 
decreasing function of r. As is seen on Figure 7, a small value of r can significantly degrade 
the upper bound of power outage performance. The more the value of r decreases, the more 
quickly the power outage performance degrades. 

Then, we evaluate the impact of circuit power consumption Pc of the sensor for separate 
and time-switching receiver architectures in Fig. 8 and Fig. 9, respectively. It is seen that, 
when the density p is small (e.g., p = 0.01), the corresponding plot is a logarithm-like function. 
Specifically, the power outage probability is very sensitive to the small value of Pc, and becomes 
less sensitive when P c is larger (e.g., above 2 x 10~ 5 W). For example, the upper bound of 
power outage probability for PPP increases from 15.1% to 58.4% (i.e., 43.3% difference) when 
Pc varies from 0.2 x 10" 5 W to 0.7 x 10 5 W. However, when Pc varies in the same amount 
from 2.0 x 10 _5 W to 2.5 x 10 _5 W, the upper bound only changes from 82.9% to 86.1% (i.e., 
3.2% difference). This indicates that, in practice, advanced sensor circuit with small Pc can 
help to lower down power outage probability significantly, especially when the available RF 
energy harvesting rate is small. For the scenarios with larger density p (e.g., p = 0.03), the 
outage probability tends to grow smoothly with the increase of Pc- The gap between the upper 
bounds of the DPP {a = — 1) and the PPP increases with the density p. This can be understood 
since a larger number of random RF sources results in a larger variance in RF energy harvesting 
rate (shown in (14) in Theorem 1), thus a greater difference in the power outage probability. 
For time-switching, as shown in Fig. 9, the impact of Pc on the power outage probability is 
dependent on the coefficient r. The larger r is, the more dynamically the upper bound of power 
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Upper Bound of Power Outage Probability versus Circuit Power Comsumption 



Upper Bound of Power Outage Probability versus P Q and i 



Fig. 8. Upper bound of power outage probability versus Fig. 9. Upper Bound of Power Outage Probability versus 
density of ambient RF energy source (separated receiver Circuit Power Consumption and Time-Switching Coefficient 
architecture). r (time-switching architecture). 


outage probability varies when P c is small. 

Lastly, we emphasize that in practice, the computation of the upper-bound obtained in Theo¬ 
rem 2 is very fast. Therefore, we believe that these bounds are the preferred choice in practical 
applications. 


C. Upper Bound of the Transmission Outage Probability 

Next, we give a practical upper bound for the estimation of the transmission outage probability 
P to . Again, we interpret these upper-bounds in terms of a worst-case scenario, which was 
specified in Section II-C3. 


Theorem 3. Let us define: 


a A / PsqPGsGh (h 0 + (1 


P c h 0 + pa 2 (2 m /bi w ) - 1) 


( 20 ) 


Then we obtain 


P to = F(C <m) < ( JJ (1 + a 


\n> 0 


T (ji T 1,7rpinf(f?, 


n\ 


-l/a 


( 21 ) 
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Proof: Applying the definition of C given in (3), it holds that 

■>m/{riW) _ -Q 


P to = P ( — [P H - ^c] + < + (2 r 

V ^ V kdK ) 


= F iP u < P c + r] (a 2 + £P H ) 


2 m/(riW) _ ^ 


= P (P H (ho - (2 m/ ^ w) - 1)) < h 0 P c + pa 2 (2 m/ ^ w) - l)) 


where we have used the fact that 2 m ^ vW ' > — 1 > 0. This implies that if /i 0 — £ ( 2 m ^ vW ) — l) < 0, 
P to = 1. Now, if ho — £ — l) > 0, it remains to reorganize the equation in order to use 

Theorem 2: 


P to = P P H < 


h 0 P c + pa 2 (: 2 m /G w ) - l) 

h 0 - rf (2 m /^w) - 1) 


( 22 ) 


so we conclude by Theorem 2 that P to < ^flnx) (-*- + ) 


1/a 

, where 


Pm 


A 

An 


\ 


gfiPsGsGu 


h 0 P c +rt(T 2 (2WO7WO-1) 
h 0 -rii[2 rn /^ w )-l) 


I PsePGsGn (h 0 + ^ (1 - 2™/(^)j) 
47 t y P c h 0 + pa 2 (2 m d r i w ) — 1) 


We note that equation (22) is an equality, so the only inequality in the proof of Theorem 3 
traces back to the application of Theorem 2. So similar to the observations made after the proof 
of Theorem 2, we note that the probability appearing on the r.h.s. of (21) corresponds to the 
worst-case scenario introduced in Section II-C3. Estimation of the upper-bound by simulation 
will therefore consist in assuming the worst-case scenario. 

We may also estimate the transmission outage probability when /C is a Poisson process as 
follows. 

Corollary 2. Let K, ~ Poiss(0,p) be a Poisson process on O = 13(0, R) with density p. Then, 
the following bound holds: 

Pto < e -^ inf ^ m ) 2 , (23) 


where p m is defined in (20). 
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Upper Bound of Transmission Outage Probability versus Density of 
Ambient RF Energy Sources (out-of-band transmission) 



(a) Out-of-band Transmission (d = 50, m = 3) 


Upper Bound of Transmission Outage Probability versus Density of 
Ambient RF Energy Sources (in-band transmission) 



(b) In-band Transmission (d = 5, m — 0.02) 


Fig. 10. Upper bound of transmission outage probability versus density of ambient RF energy sources. 


Upper Bound of Transmission Outage Probability versusp and x 



Fig. 11. Upper Bound of Transmission Outage Probability versus Density of Ambient RF Energy Sources and Time-Switching 
Coefficient r (out-of-band transmission). 


In Fig. 10, we evaluate the analytical expressions of the upper bound of transmission outage 
probability for both scenarios of out-of-band and in-band transmission. For the former and latter 
scenarios, we set the distance between the sensor and data sink d to be 50m and 5m, and the 
transmission rate requirement to be 3kbps and 0.02kbps, respectively. The numerical results are 
averaged over 10 6 runs of simulation. Similar to the upper bound of power outage probability, the 
upper bound of transmission outage probability demonstrates a similar pattern, i.e., a decreasing 
function of density p. It is shown that the analytical upper bounds for both receiver architectures 
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Upper Bound of Transmission Outage Probability versus Transmission 
Rate Requirement (out-of-band transmission) 



Upper Bound of Transmission Outage Probability versus Transmission 
Rate Requirement (in-band transmission) 



(a) Out-of-band Transmission (d = 50) 


(b) In-band Transmission (d = 5) 


Fig. 12. Upper bound of transmission outage probability versus transmission rate requirement (separated receiver architecture). 


in both scenarios of out-of-band and in-band transmission are very accurate. 

For the time-switching architecture, we examine the impact of the coefficient r in Fig. 11. We 
can observe that the larger the density p is, the more the value of r affects the transmission outage 
performance. Unlike power outage probability which is a decreasing function of r, transmission 
outage probability has a minimal value attained for r 6 (0,1). For a certain density p > 0, 
when r varies from 0 to 1, the upper bound of transmission outage probability first decreases, 
but begins to increase quickly after r exceeds a certain value. Then, after r reaches another 
certain threshold, the upper bound remains to be 1. This implies that there exists some tradeoff 
between the energy harvesting time and data transmission time to minimize transmission outage 
probability. Another interesting finding in Fig. 11 is that the optimal value of r to minimize 
the transmission outage probability are not dependent on the density p. This will be proved in 
Section III-D. 

Next, in Fig. 12 and Fig. 13, we illustrate the impact of transmission rate requirement on 
the upper bounds of transmission outage probability for separated and time-switching receiver 
architecture, respectively. It is clear from the figures that the transmission outage probability 
is an increasing function of transmission rate requirement for both architectures. As can be 
seen from Fig. 12, a larger density p causes a bigger difference between the upper bounds of 
transmission outage probability for the DPP (a = — 1) and the PPP This is because when the 
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Upper Bound of Transmission Outage Probability versus Upper Bound of Transmission Outage Probability versus 

Transmission Rate Requirement andr (out-of-band transmission) Transmission Rate Requirement andx (in-band transmission) 



(a) Out-of-band Transmission (d = 50, p = 0.03 ) (b) In-band Transmission (d = 5, p = 0.03 ) 

Fig. 13. Upper bound of transmission outage probability versus transmission rate requirement and r (time-switching 
architecture). 


density p is small, there are few RF sources, and their correlation structure does not have such 
a high impact. Similar to the study in power outage probability, the difference when a varies 
is caused by the increased variance of RF energy harvesting rate. Indeed, since the variance is 
higher, when a is close to —1, the transmission outage probability is also higher. For out-of- 
band transmission, the upper bound of the transmission outage probability is a sigmoid shape 
function of transmission rate requirement. Specifically, when the transmission rate requirement 
is relatively small or high, the upper bound of transmission outage probability increases more 
slowly with the transmission rate requirement, compared to the case when the transmission rate 
requirement is median. By contrast, for in-band transmission, the transmission outage probability 
is close to a linear function of rate requirement. In Fig. 13, for both scenarios of out-of-band 
transmission and in-band transmission, we can observe that the optimal r is dependent on the 
rate requirement. Generally, the larger the rate requirement, the smaller the value of optimal r. 
This motivates us to study the optimal value of r in the next subsection. 

Based on Theorem 3, we can also derive the lower bound of achievable transmission rate of 
the sensor. For brevity, the characterization of the lower bound of transmission rate is provided 
in Appendix II. 
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D. Optimal Choice of the Parameter in the Time-Switching Architecture 


In this subsection, we consider the choice of the parameter r in the time-switching architecture. 
Let us recall that the time-switching architecture corresponds to the choice of q = r and r/ = 1—r. 
First note that the bound of Theorem 2 is a decreasing function of r. In other words, in the 
worst-case scenario, the power outage probability is lowest when r is largest. 

Turning our attention to the result of Theorem 3, we note that the bound obtained therein is 
equal to 1 for r = 0, and also equal to 1 for r = 1. It follows that the upper bound to the 
transmission outage probability is minimized for a certain r in [0,1]. Therefore, the value t* 
which minimizes the transmission outage probability is given by 


t* = argmin re[0>1] 




^ T(n + 1, Tip mf(R, 7 m ) 2 ) ^ ^ 



(24) 


where we emphasize that in (24), 7 m depends on r. Noticing that the function 


r 


n> 


\n> 0 


a- 


T(n + l,7rpinf(f?,7 m ) 2 ) 


n\ 


-l/cx. 


is decreasing, it follows that t* can be obtained from the simpler expression 

t* = argmax re[01] 


r (ho + (1 - t)£ (1 - 2 m /(( 1 " T ) M/ ))) ) 
P c h 0 + (1 - r)cr 2 (2 m /((i-L>wO - 1) J ’ 


(25) 


where in the last expression it becomes clear that t* depends only on h 0l £,m,W, Pq and a. 
Note that, for a certain scenario when h 0 , £, W, Pq and a are fixed, the optimal r is only 
dependent on m. This analysis validates the observation in Fig. 11 that optimal r is not affected 
by density p. In addition, the analysis agrees with Fig. 13 that the rate requirement influences 
the optimal value of r. 

In Fig. 14, we present the numerical results for the upper bound of transmission outage 
probability versus r, under various rate requirement m for both scenarios of out-of-band (£ = 0) 
and in-band transmission (£ = 1). The results are averaged over 10 6 runs of simulation. It 
is obvious that the optimal r that minimizes the transmission outage probability differs based 
on m and £. This agrees with the above analysis. Then, we compare the optimal values of r 
observed from Fig. 14 with the corresponding analytical results obtained from (25) in Table 
III. The analytical results can provide very accurate guideline for the setting of optimal r to 
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Upper Bound of Transmission Outage Probability versusx 



Fig. 14. Upper Bound of Transmission Outage Probability versus Time-Switching Coefficient r (p = 0.05, a = — 1). 

TABLE III 

Optimal Values of t. 



Analytical result 

Simulation result 

£ = 0 

m = 2 

0.6828 

0.68 

m = 4 

0.4364 

0.44 

m = 6 

0.2690 

0.26 

£ = 1 

m = 0.02 

0.9185 

0.92 

m = 0.06 

0.8658 

0.86 

m = 0.08 

0.7980 

0.80 


perform best in a worst-case scenario. For discrete time-switching [46], the analytical results 
can be approximated to the closest practical value of r. 

IV. Conclusion 

This paper has presented the performance analysis of a wireless sensor powered by ambient 
RF energy through a stochastic geometry approach. We have analyzed the general cases when 
the ambient RF sources are distributed as a Ginibre a-dctcrminantal point process (DPP), which 
covers the Poisson point process case when a approaches zero. We have derived the expression 
of the expectation and variance of the RF energy harvesting rate. We have further characterized 
the worst-case performance of the sensor node using the upper bound of power outage probability 
and transmission outage probability. Additionally, we have studied the optimal value of the time- 
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switching coefficient to minimize the transmission outage probability. Numerical results validate 
all the analytical expressions, which leads us to believe that these analytical expressions are 
usable in practice. We have found that given a certain network density, the sensor achieves 
better performance when the distribution of ambient RF sources shows stronger repulsion and 
less attraction. Our system model can be extended by considering multiple-input and multiple- 
out communication channel between the sensor and data sink. Another direction of our future 
work is to extend the performance analysis from an individual node level to a network of nodes. 


Appendix I 

Proof: Recall the following fundamental formula of point process theory found, e.g., in 

[36], 

V Ph = [ (P*) 2 pU(k)dk+ [ P«p ( 2 ) (M)dMZ- ( f P*pW(k)dk) , 

Jo J OxO \Jo ) 

where p (1) and p <:2) are the first and second correlation functions defined in (5). Noting that 
p( 2 \k, l ) = K(k, k)K(l, l ) + a\K(k, l) | 2 and recombining the terms, one finds 

Var (P H ) = / (PuY K(k,k)dk + a [ P^P l n \K{k,l )\ 2 dkdl. 


<o 


’0x0 


It remains to substitute the expressions of P{\ and A'(7c, l) to obtain 


V Ph = gfJPs 


GsGrA 2 


(47t) 2 
—2 ) pdk + a 


0 x 0 (e + 11Xfc11) (e + ||x/||) 


(|fc| 2 +W 2 ) 


dkdl 


lo V(e + ||x fc | 

By polar change of variable in the first term, and simply rewriting the second, we have 


V Ph = pf3P s 


G S G U X 


2 \ 2 


(47t) 2 

2np 


»R 


dr + ap 2 


= -7rp||x-y|| 2 


Ioxo (e+ ||x||) 2 (e + ||y||) 2 


I o (e + r ) 4 

We conclude by simply evaluating the first integral analytically and obtain (14). 


dxdy . 
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Appendix II 

This appendix derives the worst-case expectation of the transmission rate of the sensor, based 
on Theorem 3. 


Theorem 4 . The expected transmission rate of the sensor node is lower-bounded as follows: 

— l/a^ 


E[ C ] > sup M 1 - TT 1 

Me(0,+oo) 1 ' ■ LA ' 


yn> 0 


a- 


T(n + l,7rpinf(f?,7 M ) z 


n\ 


(26) 


where 7 m was defined in (20) which we recall here: 


A IPsqPGsG h (ho + ^ (1 - 2 M /(nW))) 

7M 4tt y P c h 0 + pa 2 (2 M /(vW) _ 1 ) 

Proof: Since C > 0, we apply Markov’s inequality for any M > 0, 


(27) 


E [C] > MP (C > M ), 


and it suffices to use Theorem 3 to bound the probability appearing on the r.h.s.: 

— 1 /o' 


E[C] > M 1 - n 1 


\n> 0 


a- 


T(n + l,7rpinf( J R,7 M ) 2 ) 


n\ 


(28) 


Lastly, it suffices to remark that the r.h.s. of (28) tends to zero as M tends to zero and as M 
tends to infinity. Since as a function of M, the r.h.s. of (28) is continuous, the supremum over 
all M e (0, Too) is therefore finite. ■ 
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